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In this paper we derive explicit expressions for the transient state 
probabilities of the Kendall birth-death process, with and without 
immigration, for any initial condition. We then propose this process 
as a model for special services point-to-point demand, in which the 
births represent circuit "connects" and the deaths represent "discon- 
nects. " This choice of model is based on intuitive arguments and on 
the fact that the model can represent the growth and turnover 
characteristics of special services demand. Thus, the model provides 
a means by which special services demand, with its inherent uncer- 
tainty, may be approximately represented in various facility network 
studies, to obtain, at the very least, useful qualitative results. In 
particular, we evaluate the probability of a held order (i.e., the 
probability that a service request is held for lack of spare facilities) 
with Blocked Customers Held (bch) as the queue discipline. We also 
apply the model to capacity expansion problems, introduce the con- 
cept of margin, the extra capacity needed to meet the demand within 
a given held-order probability, and examine its sensitivity with re- 
spect to growth, turnover (or churn), and system size. We find that 
aggregating small demands into a single larger demand produces 
significant reduction of the margin, because of improved statistical 
properties. 

I. INTRODUCTION 

In this article, the transient behavior of the Kendall birth-death 
process* with immigration is examined, and some applications of the 



* The Kendall birth-death process is one in which the transition rates are proportional 
to the state. 
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process to capacity expansion problems are discussed. The choice of 
such a process was motivated by the search for a model for special 
services point-to-point circuit demand, a model which would be used 
as a tool for determining facility network circuit routing strategies. 
Special services demand generally consists of demand for full-time 
dedicated circuits (e.g., foreign exchange lines, wats lines, data lines), 
as opposed to the message-traffic offered load which consists of de- 
mand for the use of common facilities for a relatively short period of 
time. Thus, the system examined is characterized by the stochastic 
process Mf) with realizations (states) n = 0, 1, •••,00, where n might 
refer to the number of working circuits or some other facility, rather 
than to the number of busy trunks, as in the message-traffic case. By 
definition, the birth-death process 1 allows transitions from some state 
n to n + 1 via a birth (circuit connect), or to n — 1 via a death (circuit 
disconnnect). The transition rates are A„ for the births and /x„ for the 
deaths, both of which are chosen proportional to n for the following 
reasons. 

It is clear, for special services, that the rate of disconnects, /*„, is 
state dependent. There are, in fact, indications 2 that ju„ is a monoton- 
ically increasing function of n. The simplest such function is n/x, which 
implies that the probability of disconnects is proportional to the size 
of the system. With this choice for the death rate, a number of possible 
choices exist for the birth rate. Choosing it to be a constant causes the 
mean number of circuits to saturate in time, while choosing it to be 
proportional to n causes the mean to grow or decline exponentially. 
Since special services are presently characterized by significant net 
growth, it would seem that a plausible model for special services 
demand is a birth and death process in which both the birth and death 
rates are proportional to the state. 

One consequence, however, of assuming A„ = rik is that if the process 
reaches the state n = at any time, by a succession of disconnects, it 
will stay there forever, since the birth rate is zero. To eliminate this 
characteristic, the concept of immigration may be introduced by taking 
\ n = n\ + (S, where /? is the immigration factor. The cases with and 
without immigration will be discussed below. 

It must be emphasized that it is not the intent of this paper to 
validate the model based on an examination of actual special services 
data. Such statistical data analysis is important for a final assessment 
of the accuracy of the model and is currently being undertaken. For 
the purposes of this paper, it shall be assumed that a study of the 
proposed model is justified, based on the intuitive arguments given 
above and on the fact that the model captures the growth and turnover 
characteristics of special services (see Section 5.1). The model provides 
a means by which special services demand, with its inherent uncer- 
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tainty, may be approximately represented in various special services 
facility network studies, to obtain, at the very least, useful qualitative 
results. 

Since the situation of interest is that of net growth, it is clear that 
statistical equilibrium does not exist, and that it is the problem of the 
transient solutions of the Kolmogorov birth-death equations that is of 
prime importance. Much literature exists on the subject of transient 
solutions for birth-and-death processes 1 ' 3 " 11 and the case in which the 
transition rates are state independent is completely solved. 7 ' 8 The case 
in which A„ and n„ are proportional to the population is solved when 
immigration is not included: The results for a specific initial condition, 
namely starting from the state n = 1, are derived in Refs. 4 and 10 and 
the expressions for the general initial condition are quoted in Ref. 10. 
For the nonzero immigration case, the form of the generating function 
for the state probabilities is known, 10 but it seems that explicit expres- 
sions for the state probabilities have not previously appeared in the 
literature. In this paper, these expressions are derived for any non- 
negative value of /?. 

In special services, if an order for service is delayed because of lack 
of spare facilities, the order is said to be held. Thus, in order to study 
capacity expansion problems, the probability of a held order is intro- 
duced, as well as the concept of margin, the extra capacity needed to 
meet the demand within a given held-order probability. This held- 
order probability is similar but not identical to the transient time 
congestion of the process (see Appendix B). The queue discipline 
followed here is Blocked Customers Held (bch), in which an arriving 
customer spends a total time T (random variable) in the system, after 
which he departs regardless of whether he is waiting to be served (i.e., 
his service order has been delayed) or is actually being served (i.e., he 
has been assigned a circuit). 

A fundamental difference between this analysis and teletraffic must 
be emphasized. This difference arises because of the respective time 
scales in the two cases. Whereas the mean lifetime of a call in traffic 
(t = l/n) is of the order of a few minutes, the mean lifetime of a circuit 
in the process described here is of the order of a few years. It is this 
fact, coupled with the relatively fast growth of special services demand, 
that makes it impossible to even approximately treat the process in a 
statistical equilibrium mode (no growth) with a slowly varying enve- 
lope representing the growth. Thus, the transient aspect of the problem 
is to be contrasted to the more conventional assumption, in teletraffic 
theory, that statistical equilibrium prevails (it must be mentioned, 
however, that some work has been done concerning nonstationary 
telephone traffic with time-varying Poisson-offered load, e.g., Refs. 12 
to 14). It must be further noted that, although the model is being 
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proposed for special services demand, nevertheless, it may be applied, 
with an appropriate choice of parameters A, /?, and /i, to any process 
that behaves in a similar manner. 

The held-order probability having been defined and the concept of 
margin introduced, questions concerning capacity expansion problems 
are addressed. Capacity expansion is a problem that has been studied 
by many. In this paper, optimal capacity-expansion policies are not 
sought; only very specialized aspects of the problem are considered. 
For instance, the effects of aggregating demands into a larger single 
demand are examined, and the minimum capacity increment which 
would meet the demand within a specified interval of time and within 
a given held-order probability is determined. In addition, the relation- 
ship between spare capacity and lead time is discussed (see summary 
of results in Section II). Some relevant work has been done by 
Freidenfelds 1516 in which the author computes first-passage times to 
various levels of demand using a general birth-death process, and 
discusses briefly fill-at-relief problems. Work by Luss and Whitt 17 
studies the impact of both deterministic and stochastic models on 
utilization. The authors use Brownian motion to model the stochastic 
demand and follow a scheme similar to ours for deterniining the margin 
needed at a future time. 

The organization of this paper is as follows. Section II sets up the 
problem and gives a summary of results. The explicit solutions for the 
general case are derived in Section III, and their properties are exam- 
ined in Section IV. In Section V, growth, turnover, and churn are 
defined, the concept of margin is introduced, and some of its applica- 
tions to capacity expansion problems are discussed. Finally, Section 
VI contains the conclusions. 

II. BACKGROUND AND SUMMARY OF RESULTS 
2. 1 General birth-death equations 

Consider a system described by a set of states n = 0, 1, •••,<», and 
a birth and death process defined by a set of transition rates {\„, ju„} . 
The quantity A„5(ju„5) + o(8)* is the probability of a birth (death) in 
the small interval [t, t + 8], given that the system is in state n at time 
t. 1 The probability of more than one birth or death in [t, t + 8] is o(8). 
The probabilities P n (t) of finding the system in state n at time t must 
satisfy the well-known infinite set of difference-differential equations 
(p. 454ofRef. 1) 

d 

-J-Pnit) = -(Xn+ Hn)Pn(t) + \ n -lPn-l(t) + jU„ +1 P„ +1 (f) 

at 

[for n>0, P-i(*)=0, mo = 0]. (i) 



• o(-): J? 1 -» /?' is such that lim ^ = 0. 
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If the initial number of circuits is no, the initial condition may be 
written 

Pn(0)= Sn^ (2) 

where S n/l0 is the Kronecker delta. 

The particular birth-death processes considered in this paper are 
the cases in which the transition rates are proportional to the popu- 
lation, n, with or without immigration. 5,610 The corresponding transi- 
tion rates, defined for all nonnegative integers, n, are 

X n = nX + /?, [in = nil, (3) 

where X, ju, and (5 are nonnegative constants. In the following sections, 
results for the case with no immigration may be easily obtained by 
setting P = 0. 

2.2 Mean and variance 

It has been shown 1,410 that the mean, m(t), and the variance, v(t), of 
processes such as those described by eqs. (1) and (3) may be obtained 
without solving explicitly for the P n (t). The resulting expressions, 
satisfying initial condition (2), may be easily found to be 

(i) Case X ^ /t 

m(t) - ( no + t-£-V~* ,< " T^-> < 4) 

\ X-fiJ X-n 

X-n 

+ - P „ [\e**-> )l - {X + n)eW + /i]. (5) 

(X - ilY 

(ii) Case X = y. 



m(t) = (St + no, 


(6) 


v(t) = X(S? + (ZXno + P)t. 


(7) 


2.3 Solutions for P„(t) 




To simplify the notation, define the following quantities: 
A = X - /x, 




AV 




A () iiMe" - l) v 




e A '-l 
R R(t\ 




Aer — p 




A 


(K\ 



w ww Xe*' - m' 
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and 



r(r-l) • • • (r - m +1) 
m! 



(9) 



This definition of the binomial coefficient is valid for any real number 
r and any positive integer m (see p. 50 of Ref. 1). For m = 0, one defines 
(S) = 1, and for negative integers m, one defines (m) = 0. The symbol 
(m) is not used if m is not an integer. Denoting v = /?/A, where v is any 
nonnegative real number, the solutions derived in this paper are 
(i) Case \¥^ fi 



Pn(t) = 



C'(/iB) no (AB)" £ 

t-0 



/lo + n + p — I — 1 

n — i 



(A - 1)', if no + v > 0, 
[S„o, if «o + v — 0. 
(«) Case A = ju 



(10) 
(ID 



P»(*)- 



A* 



1 + \t) VI + A* 



n+ "o min</io,/i) 

2 

i-0 



r) 



rio + n + p— i— 1 

n — i 

8 n0 , if no + v = 0. 



^-1 I, if no + v>0, (12) 

(13) 



Equations (10) and (11) with no immigration (v = 0) are identical to 
the results quoted by Bailey [Eqs. (8.47) of Ref. 10]. 



2.4 Application to capacity expansion 

In Section V, margin is defined as the capacity which must be built 
in excess of the mean to meet certain service requirements, and the 
percent margin is defined as the ratio of the margin to the mean in 
percent. The following is a summary of the main results: 

(i) By aggregating demands, less percent margin is needed than in 
the nonaggregated case. This effect is especially significant for small 
demands. 

(ii) Given a niinimum desired time, T, between successive expan- 
sions, a procedure is established for determining the minimum capacity 
increment which would meet the given service requirements. 

(Hi) Given a lead time, t, between the moment faculties are ordered 
and the time they are available for use, a procedure is established for 
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determining the threshold value of the remaining spare corresponding 
to the time at which new facilities should be ordered. 

(iv) By introducing immigration, the absorbing zero state is elimi- 
nated and the percent margin needed to meet the service requirements 
is reduced for moderately large to large times (of the order of two years 
or more for the particular values examined). 

III. DERIVATION OF THE STATE PROBABILITIES 

The approach followed to solve the set of equations in (1) is the 
generating function technique. 510 In Ref. 10, a differential equation for 
the generating function, F(s, t), defined below, is established and its 
solution is derived. The results are quoted in Section 3.1. Three well- 
known identities are given in Section 3.2 and are then used in Section 
3.3 to derive explicit expressions for the state probabilities. The pro- 
cedure followed in Section 3.3 is to identify F(s, t) as the generating 
function for a convolution of two functions. 

3. 1 The generating function 

The generating function, F{s, t), is related to the state probabilities 
through the following expression: 

F(s,t) = l s n P n (t), 0<s<l. (14) 

77=0 

The differential equation for F(s, t), given in eq. (8.63) of Ref. 10 with 

e e = s, is 

«£« + ,„.,«£«>.,<._ ,,,<*,>, (1S ) 

dt ds 

where 

H(s) = -(s-l)(As-/x). 
The solutions to this equation are 

A Y/b + asY 



F(s, t) = 



d + csj \d + cs 
b + as 



d + csj \d + cs 



A + li, (16) 

A = /C (17) 



* An alternative approach for obtaining this result is to substitute \ — A for n in the 
A^n expression and to take the limit A — » 0. 
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where 

< s < 1, 

a = A - ne* 1 , 

b - -ii(l - e"), 

c = A(l - e A '), 

d = Ae At - n, (18) 

d = 1 - \t, 

b = \t, 

c = -\t, 

d=l + \t. (19) 

The above solutions may be verified by direct substitution. Equation 
(16) agrees with eq. (8.71) of Ref. 10 and eq. (17) with v = agrees with 
eq. (8.52) of Ref. 10. 

3.2 Useful identities 

In Section 3.3, use will be made of the three following well-known 
identities. 

3.2.1 Binomial identity 

For any a and /? and for any nonnegative integer n, the following 
identity holds (p. 51, Ref. 1): 

(« + £)"= i [ " |a"-T. (20) 

m -0 \ m / 

3.2.2 Negative binomial identity 

For any a and ft such that | fi/a \ < 1 and for any real number r, the 
following identity holds (see pp. 51 and 269 of Ref. 1): 

(a-pr r = V (-l) m ( r )B m a- {m+r) . 

m-o \ m J 

If r is strictly positive, identity (12.4) on p. 63 of Ref. 1 may be used to 
write 

r + m — 1 \ nm _,„, + r! 



m=0 



(a-i8)- r = SI )0 m a- lm+r \ (21) 



m 



3.2.3 Generating function for a convolution 

Let /^(s) and ^2(5) be the generating functions for the sequences 
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{P^'Jn-o.i,- ■ • and {Pf }„-o.i,-- •. respectively, 

P,(«) = £ s*P£\ F 2 (s) = £ s"Pj?\ (22) 

The function F(s) = Fi(s)F 2 (s) is then the generating function for 
{Pn} n -o,i," •. the convolution of Pi" and Pn\ and may be written as 

F(s) = l s n P n , (23) 

71-0 

where 

Pn = l PFP& = S PPW2.-. 

«=0 i=0 

The proof of this theorem is elementary (e.g., see Chapter 11 of 
Ref. 1). 

Note: This theorem applies to arbitrary sequences {Pn ] } and {P™} 
(not necessarily probability distributions) as long as their respective 
generating functions exist. Thus, the series in eqs. (22) must converge. 
For the purpose of this theorem, however, it is assumed that Fi (s) and 
F2(s) do exist. 

3.3 Derivation of explicit expressions 
3.3.1 Case A ^ /t 
The generating function in eq. (16) may be rewritten as follows: 

F(8, t) - Fi(s, t)F 2 (s, t), (24) 

where 

Fi(s, t)=A"(d + cs)- lno+v \ 

F 2 (s, t) = (b + os) n °. 

(a) rio + v>0 

Applying identity (20), it may be seen that F 2 (s, t) is the generating 
function for a binomial type function, 

F 2 (s,t)= I ("° )b n °- m (as) m 
m -o \tn I 



= £ s-P<?(«, (25) 

m=0 
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where 



(Z)/*\ = . 



P%(t) = 



if m < no, 



if m > no. 



(26) 



In a similar manner, it may be seen from identity (21) that F\(s, t) is 
the generating function for a negative binomial type function, 

F 1 (s, O-A'2 (*o + " + >*- ^(-cs)^-""— ' 

m=0 \ m / 



where 



= £ s-Pii'U), 

m = 



P ( t ) = A*( Bo + " + m _ ^ ^ " ' '" ' ; : "'' ' 



(-c) m d 



(27) 



(28) 



It may be shown that | cs/d | < 1 for all values of t > 0, < s < 1, and 
A, n > 0. Thus, identity (21) applies in all the relevant cases. 

It now follows from eq. (23) that F(s, t) is the generating function of 
the convolution, 



Pn(t) = I P\»(t)P™i(t) 

j-0 

min {/i () ,n ) 
= I 



no 



no + v + n — i — 1 



/i - * 



(-c)"-'d 



n-i J — (n a +v+n — i) 



, (29) 



where the upper limit on the sum arises from the condition i < no for 
Pi 2) (t) established by eq. (26). 
Rearranging, one obtains 

(A \ " min (n .n) 
3 ,? 



no\(no + n + v — i — 1\( a \ ( & \ 



n — i 



(30) 



From the definitions of a, b, c, and d in eqs. (18) and from eqs. (8), the 
above expression for P n (t) reduces immediately to eq. (10). 

(b) no + v = 

For this case, F(s, t) = 1. From the definition in eq. (14), it is then 
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apparent that all P n (t) for n ¥^ must be zero and that P (t) - 1, 
which is the result shown in eq. (11). 

3.3.2 Case A = /t 

The solutions may be obtained by using the same procedure followed 
in Section 3.3.1. The only difference is that the starting equation 
should be eq. (17) rather than (16). Since eq. (17) can be simply 
obtained from eq. (16) by letting A — > 1, a -* d, b — > b, c — > c, and d 
— > d, it follows that the final results for the case A = /i can be obtained 
from the results of the case A # fi [i.e., eq. (30)] by making the above 
substitutions. Alternatively, eqs. (12) and (13) may be obtained by 
taking the limits of eqs. (10) and (11) as /z — » A. The procedure is the 
following: Consider A to be a constant, then replace /i by A — A wherever 
it appears, and finally take the limits A — > using l'Hospital's rule 
whenever necessary. The results of this limiting procedure are found 
to be 

/<— A \At 

hmB{t)= T^T7' < 31 > 

M — A 1 + At 

limC(*) = 



^x 1 + Xt 



IV. PROPERTIES 



In this section, the zero-state probability and the cumulative prob- 
ability distributions, for several choices of the parameters A, ju, and fi, 
are examined as a function of time. In addition, some cases in which 
the state probabilities are especially simple are indicated. 

4. 1 Probability of ultimate extinction 

The birth-death process without immigration is characterized by an 
absorbing state at n = 0. If the system reaches that state at some time 
to, it will stay there for all t > to since the birth rate is zero. The 
probability of hitting that state at time t is given by Po(t). In the limit 
t — > oo, this probability tends to 



limP (0 = 

t—oo 



gl. *>, 



(32) 
1, A < /i. 

Thus, for any fi^O, there is a nonzero probability of ultimate extinction 
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Fig. 1 — Zero-state probability. 



of the process, while for ju > A, ultimate extinction is a certainty. This 
feature may be removed, if desired, by introducing immigration. In 
this case, the limiting value of Po(t) tends to a nonzero value for A < 
jtx and to zero for A > /i, 



lim Pott) = 

/-.00 



x) ft)' £=•""-«• x> * 



i-Jl. 



(33) 



The effect of immigration on Pott) is shown graphically in Fig. 1, 
where Pott) is plotted for various values of A, ft, and v, and for no = 1. 
The case no = 1 was chosen for clarity of the figure, since the effect of 
immigration is larger for smaller values of no. 

4.2 Cumulative probability distribution 

The cumulative probability distribution is defined as 

Pntt) - l Pdt). 

i-O 

In order that \im n -,coF n (t) = 1 for all t, it is necessary and sufficient 
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that X*=o A^ 1 diverges, which is the case for the process discussed in 
this paper (see Theorem on p. 452 of Ref. 1). The function, F n (t), for 
various choices of X, ju, and ft, and for no = 5 is shown in Figs. 2 through 
5. The lowest curve plotted in each figure is Po(t), and is consequently 
a measure of the extinction probability. 

The values of A (0.6) and /x (0.3) chosen in Figs. 2 and 3 correspond 
to net positive growth (see discussion of growth in Section 5.1). As 
may be verified from eqs. (32) and (33), the lim,_»P U) is nonzero in 
Fig. 2 and zero in Fig. 3. In addition, for each n > 0, the lim^=cPnU) 
= 0, although the lim,^ £"-i P n (t) is nonzero. Thus, for any n > 0, 
the lim,^ooF„U) is nonzero for ft = and zero for ft ^ 0, reflecting the 
fact that the extinction probability is nonzero in the first case and zero 
in the second. (The n = curve in Fig. 3 is essentially flat and is hard 
to distinguish on the graph.) 

The case in which the death rate is larger than the birth rate is 
shown in Figs. 4 and 5. This situation corresponds to negative growth. 
As may be verified from eqs. (32) and (33), the lim,_oc.PoU) is 1.0 in 
Fig. 4 and 0.5 in Fig. 5. It may be shown that, for the ft = case, 
lim,_ocP„(0 - for all n > 0, while for the ft ^ case, lim,^«P„U) * 
for all n > 0. In both cases, however, the lim,_ooF„(0 is nonzero, 
reflecting the fact that the extinction probability is nonzero. 

Finally, the case A = /n, ft = 0, is similar to Fig. 4 (with extinction 
probability equal to unity), and the case A = ju, ft 5* 0, is similar to Fig. 
3 (with extinction probability zero). These cases are not shown. 
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Fig. 2 — Cumulative distribution without immigration (A > n). 
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Fig. 3 — Cumulative distribution with immigration (X> fi). 

4.3 Special case solutions 

For certain initial conditions, the general solutions reduce to simple 
analytical forms. For no = 0, the interesting process is the one including 
immigration (p > 1). The state probabilities of eq. (10) may then be 
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Fig. 4 — Cumulative distribution without immigration (a < u). 
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Fig. 5 — Cumulative distribution with immigration (a < fi). 



written 



PM) = C'(\B) n 



n + v- 1 

v- 1 



(34) 



For no = 1, the process without immigration becomes interesting. 
Equation (10) with v = reduces to the well-known form 410 

P (t) = fiB, 

P n (t) = biB)(\B) n A = (1 - AB)(1 - /iflMAB)"- 1 



(n ¥> 0). (35) 



V. APPLICATIONS TO CAPACITY EXPANSION 

In this section, the A, /x, and /? parameters of the model are related 
to more physically intuitive quantities such as growth and turnover 
(or churn). The concept of margin, the extra capacity needed to meet 
the demand within a given held-order probability, is introduced. The 
effects of randomness and immigration on the margin are then exam- 
ined, and finally, several capacity expansion problems are addressed. 

5. 1 Growth, churn, and turnover 

The model described in the preceding sections is completely speci- 
fied once the parameters A, /x, and /? are known. In this section, 
quantities that are more physically intuitive than the birth and death 
rates, namely growth and turnover (or churn), are introduced and 
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related to A, ju, and /?. First define the following quantities: 

b(t) = EB(t) = mean number of births in [0, t], (36) 

d(t) = ED(t) = mean number of deaths in [0, *]. (37) 

Then, 

m(t) - no = b(t) -d(t) 

= mean net population increase in [0, t], (38) 

where E refers to the expected value and where m(t) is the mean value 
of the population, as defined in eqs. (4) and (6). Differential equations 
for b(t) and d(t) are derived in Appendix A and exact analytical 
solutions for these equations are found. 

The annual rate of growth, g, is defined as the change in the mean 
number of circuits in one year divided by its value at the beginning of 
the year. Thus, 

m(,+ l)-mW m 

m(t) 

From eqs. (4) and (6) it may be seen that 

[no + i S/(A- J u)]e (A -^(e A -"-l) 



g(t)~ 






no + fit 



for A 5^ n, 

(40) 
for A = ft. 



By observation, it may be noted that if p = 0, the growth is time 
independent, whereas if ft > 0, the growth depends on time. A time 
average value of the growth may be defined to be 

l=limi g(t)dt. (41) 

T — T Jo 

It may be easily found that 

[0, for \<fi, 0*0, 

g - 1 0, for A - ft, (42) 

Le x_ ' 1 — 1 for all other cases. 

When A = ju, it must be borne in mind that g = does not necessarily 
mean no growth. In fact, for the fi > case there is linear growth at 
the rate fi [recall eq. (6)]. Thus, it is the immigration factor that 
represents the growth in this case. 

In summary, throughout this paper, the last expression in (42) will 
be used as the definition of growth with the provision that it is the 
variable ft that actually describes the growth in the A = y. case. The 
case A < jU and /? > will be disregarded. 
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Churn may be defined in various ways. Its purpose is to quantify the 
"activity" of the process, i.e., to compare the number of "connects" 
with the number of "disconnects." For instance, churn may be defined 
as the ratio of the mean number of total connects in one year to the 
mean number of net connects in the same period. This ratio has also 
been called in-to-net, or in-to-gain ratio. 15 Thus, 

clt) ~ Ht +1) ~ W) ~ IN (43) 

K m(t+l)- m(t) net { 

This quantity has the easy interpretation of being the expected number 
of connects in one year for a net increase of one circuit. For instance, 
a churn of four (which seems to be a typical number 18 ) means that 
four connects are expected for every net increase of one circuit. Of 
course, it follows that three disconnects are also expected for consist- 
ency. The problem with this definition of churn is that for low-growth 
cases (i.e., when the net increase is almost zero) the ratio of total 
connects to net may become a very large number. Furthermore, in the 
case of negative growth, this ratio becomes negative. Thus, the range 
of values which the churn may take is very large, which makes it a 
difficult number to work with in data analysis. 

For this reason, an alternative definition of churn is introduced and 
is called turnover. Turnover is the expected number of connects 
(disconnects) needed to replace the number of disconnects (connects) 
that occurred in one year, divided by the expected number of circuits 
in place at the beginning of the year. The words outside the parentheses 
refer to the positive-growth case in which there are more expected 
connects than disconnects, and the words in the parentheses refer to 
the negative-growth case when the reverse is true. The turnover may 
be written as 

1 in + out - I in - out I 

a(t) = ■ [ 44 

2 MEAN 



min(iN, out) 



(45) 



MEAN 

where 

in = b(t + 1) - b(t), 
out = d(t + 1) - d(t), (46) 

MEAN = m(t). 

Thus, a turnover of 0.3 with positive growth indicates that over the 
next year the expected number of disconnects will be equal to 30 
percent of the mean at the beginning of the year. The expected number 
of connects depends on the growth and will be greater than or equal to 
the disconnects. (Note that it is not necessarily the same circuits that 
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are connected and disconnected.) From the expressions in Appendix 
A, eq. (44) may be written as 



a(t) = 



where 



i (/i eff e A ' - /5/A)- l [/ieirc A V - 1) (A + M)M 

- 2fli/A - n eff e A ' | (e A - 1) | ] (X * /t), 



1 2Ayfi< + A/? + 2A/io 

2 " no + /?* 



(47) 



(A - M). 



n e ff = no + 



/8 



X-/i 



Again, by observation, it may be noted that if fi = 0, the turnover is 
time independent, whereas if fi > 0, it is time dependent. As in the 
case of the growth, a time-average value of the turnover may be 
defined to be 



a = lim 



if- 



It may be easily found that 



a = 



A 
A 



M 



(t) dt. 



if X>/i, 

if \<n, /3 = 0, 

for all other cases. 



(48) 



(49) 



Throughout this paper, the above expressions for the turnover will be 
used. As mentioned before, the case A < /t and fi > will be disregarded. 

From preliminary data analysis, it has been found that typical values 
of growth and turnover for special services fall into the range —0.15 to 
+0.15 for the growth and to 1 for the turnover. Nominal values of g 
= 0.1 and a = 0.3 were chosen in this paper. 

To do a study using stochastic special services demand, numerical 
values of the parameters of the model are needed. Accurate values, if 
such values exist, may only be found by careful data analysis of 
historical demands. Approximate values, however, may be found as 
follows. First equate the mean [eqs. (4) or (6)] to the special services 
forecast to obtain values for A and fi. Then equate the variance [eqs. 
(5) or (7)] to some measure of the forecast uncertainty to determine 
A and ju. For some cases, given below, one may conveniently use eqs. 
(42) and (49) to write the mean as a function of growth alone, and the 
variance as a function of both growth and turnover. The results are as 
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follows. 

(i) Exponential growth with no immigration (A ¥* n, fi = 0) 

m(t) = n (l + g) 1 , (50) 

v(t) =n (^± 1 J (1 + g)%(l + gV - 1], 

where the positive (negative) sign refers to A > /t (A < /i), 
(«) Linear growth (A = /i, /? > 0) 

m(f) = # + n , (51) 

u(*) = a/8* 2 + (2dho + 0)*. 

5.2 Margin and minimum capacity increments 

Define the quantity h{t) as the probability of a held order, i.e., the 
probability that at least one service order is delayed due to lack of 
spare facilities. Then 

h(t)= I P n (t) = l-F d (t), 



-d+l 



where d = d(t) is the total number of servers (facilities) at time t. The 
quantity h(t) is similar but not identical to the transient time conges- 
tion function (see Appendix B). Computationally, d(t) may be deter- 
mined from the state probabilities by requiring h(t) to be less than or 
equal to some predetermined number, h. Thus 

d(*) = min{d = 0,l, ...| £ P n (t) < h}. (52) 

n-d+l 

Of course, the actual sum involved in the computation is not infinite, 
since the condition in eq. (52) is equivalent to£n-o P n (t) > 1 - h. The 
level d(t) can be viewed as the sum of the mean number of circuits 
and a quantity which may be called margin. Thus, given any time t > 
0, the margin is the capacity which must be built at to = in excess of 
the mean m(t), in order to meet the demand, within the maximum 
held-order probability, h. 

A significant quantity is the ratio of the margin to the mean in 
percent which will hereafter be referred to as the percent margin. 
Figure 6 shows a plot of this ratio as a function of time for various 
values of growth and no immigration. Turnover has been taken to be 
0.3, the initial number of circuits 5, and the maximum probability of a 
held order 0.05. As may be seen, the percent margin increases with 
time, and generally less percent margin is needed for larger growth 
rates. The sensitivity of the percent margin with respect to turnover, 
for a growth of 0.1, is shown in Fig. 7. It may be seen that the larger 
the turnover, i.e., the larger the "activity" in the network, the more 
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Fig. 6 — Sensitivity of percent margin to growth (no immigration). 



margin one has to build to provide the same maximum held-order 
probability. 

The concept of margin has several applications, one of which is the 
determination of an appropriate capacity increment at each expansion. 
Given a minimum desired time, T, between expansions, it would be 
useful to determine the minimum capacity increment, c, which, if 
installed at time t, will exhaust in [t, t + T] with a probability that is 
no larger than h, or equivalently, the increment which will last the 
interval T with a probability greater than or equal to 1 — h. Thus, the 
condition on c may be written as 

Piob{^(t + &)< no + c\^(t)** no, V £e [0, T]} > 1 - h. (53) 

Since the X and /z coefficients are time independent, the process is time 
homogeneous. Consequently, changing the origin of time does not 
affect the problem. Choosing it to be at t is equivalent to setting t = 
in the above expression, and determination of c reduces to finding 



mim c = 0, 1, 



rto+c 

I Pn(Z) 21-h, v £ e [0, T) 

n=0 



) 



(54) 
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Tables may be set up permitting direct reading of the values of c 
corresponding to the growth, turnover, initial state parameters, and to 
the time interval T. Some typical results are plotted in Fig. 8. 

For completeness, it must be mentioned that in problems of the type 
described above, questions about service-order queue disciplines must 
be entertained. A careful consideration of this aspect of the problem is 
beyond the scope of the present analysis, and the queue discipline 
implicitly followed has been Blocked Customers Held (bch). See 
Appendix B. 

5.3 Effect of randomness: Aggregation benefits 

An inspection of Fig. 6 shows that the percent margin needed is 
large, for the particular case examined. Since the initial state consid- 
ered is rather small (no = 5), an interesting question is to find out 
whether the percent margin can be reduced by combining demand to 
form larger quantities, in the hope that the statistics of the aggregated 




Fig. 7 — Sensitivity of percent margin to turnover (no immigration). 
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process will be better behaved, i.e., less susceptible to fluctuations. Our 
studies have shown that benefits are indeed obtained by aggregation. 
Figure 9 illustrates the results. The time evolution of the percent 
margin is plotted as the initial state of the system is varied, for fixed 
values of growth, turnover, held-order probability, and for no immigra- 
tion. Two important observations may be made. The first one is the 
fact that the percent margin decreases as the initial state of the system 
increases. An implication of this behavior, for example, is the following: 
Suppose demand between two points is being satisfied by two inde- 
pendent routes (with initial number of circuits no 1 ' and n 2) , respec- 
tively). Benefits would be obtained by combining the two demands on 
one route (with initial number of circuits no 1 ' + n 2) ) because the 
margin one would have to build in this case is less than in the 
nonaggregated case, the held-order probability being the same. The 
second observation is that the change in the percent margin with 
respect to n is larger for small values of no. The implication is that the 
benefits will be especially significant when aggregating small demands. 
It must be mentioned that the conclusions about aggregation 
benefits in the example given above were based on an examination of 
Fig. 9. The implicit assumption was that the combined process would 
obey the same birth and death equations as each single process, and 
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that it would be described by the same pair of transition rates, A and 
li. This assumption is only true if all the individual processes are of the 
same type, i.e., if they all have identical transition rates. It turns out, 
however, that the method for determining margin described in the 
previous sections may be extended with little additional effort to the 
general case in which there are M simultaneous processes characterized 
by a set of transition rates {A', /i'} for i = 1, • • • , M. Since the processes 
are independent, the joint probability distribution is the product of 
the individual distributions, 



Pn„n. 2 ,...,nJt)=P { n \ ) (t)P%(t) ■••Pfj(t). 



(55) 



The probability of being in some level n, regardless of the composition 
of that level, may then be written 



Pn(t) = I' I' 
n, n 2 



Z'p£!(t)P%(t)-..p% ) it), 



(56) 



n,\ + n,2 + 



+ jxm = n. 



The sums are over all values of tii, n 2 , • • •, Um such that n\ + nz + 
• • • + Km = n. The primes over the summations are an indication of 
this restriction. The margin for the combined process may then be 
determined from the mean m(t), and the quantity d(t), analogous to 
that defined in Section 5.2. The mean for the combined process is 
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simply the sum of the individual means, 
m(t) =mi(t) + m 2 (t) + 



+ m M (t), 



(57) 



and d(t) may be obtained from an expression similar to eq. (52), 
namely, 



d(t) = minj d = 0, 1, • • • 






(58) 



5.4 Effect of immigration 

Immigration affects the problem in several ways. First, it eliminates 
the absorbing state at n = 0, and consequently the probability of 
extinction, for all cases except the A < ju, /? =£ case. Furthermore, the 
lim*_oo P n (t) = (n > 0), for all cases except the case mentioned above, 
for which the limit is nonzero. Finally, a nonzero value of ft gives the 
model flexibility to represent linear growth (for A = ju) as well as 
exponential growth (for A ^ n). It is interesting to note that for the A 
> ju case, immigration actually reduces the percent margin for moder- 
ately large to large times, as seen in Fig. 10. This behavior is due to the 
fact that introducing immigration increases the mean (which tends to 
decrease the percent margin) faster than it increases the variance 
(which tends to increase the percent margin). 
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Fig. 10 — Effect of immigration on percent margin. 
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5.5 Other applications: Spare threshold 

In capacity expansion problems, a question that arises is when to 
expand. If facilities could be installed instantaneously, the expansion 
time would be simply the time at which remaining spare exhausted. 
However, there usually is a lead time, t, between the moment facilities 
are ordered and the time they are actually installed. If existing spare 
can be monitored, it would be useful to determine a priori what the 
particular level of spare would be at the time when new facilities 
should be ordered. This information would then yield the order time, 
since as soon as that spare level is attained, it is time to order. This 
value of spare, or spare threshold, a, is the amount of remaining 
capacity which will exhaust in [t, t + t] with a probability that is no 
larger than h, or equivalently, the amount of capacity which will last 
the interval t with a probability greater than or equal to 1 — h. If the 
last expansion occurred at to, providing a total capacity of C(to) (see 
Fig. 11), the constraint on a may be written as 

Prob{J^(t + £) < C(t<>) | 

JT{t) = C(to) -a, V £ E [0, t]} > 1 - h. (59) 
Time homogeneity of the process allows setting t = in the above 
expression, as discussed in Section 5.2, and the determination of o 
reduces to rinding 

mini a -0,1, ••• E P„(£) > 1 - h, V|e[0,r]L (60) 
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where a appears implicitly in the initial condition used to evaluate the 
state probabilities P n (£). Again, tables may be set up permitting the 
direct reading of the values of a corresponding to the growth, turnover, 
capacity levels, and the lead time, t. Some typical results are plotted 
in Fig. 12, where a is shown as a percentage of the total capacity. 

VI. CONCLUSION 

A summary of the main results has been given in Section II. Explicit 
solutions for the birth-death process in which the births and deaths 
are proportional to the state have been derived, and some of their 
applications to capacity expansion problems have been discussed. A 
method for determining margin has been described in detail for the 
case in which one birth-death process exists with exponential growth 
characterized by a pair of transition rates, X and ft, and it was shown 
how to extend the method to cases in which there were several 
simultaneous processes. 

In addition to its applications to capacity expansion problems, the 
proposed model is useful in assessing the potential of routing strategies 
for special services. It has already been shown, for example, that 
benefits are to be expected by aggregating small demands. These 
conclusions were based largely on service robustness considerations. 
To obtain more comprehensive results about routing strategies, it is 
clear that cost robustness considerations must be addressed as well. 
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APPENDIX A 

Expressions for b(t) and d(t) (Section 5.1) 

As indicated in Section 2.1, 

X n 8 + 0(8)* = Probability of a birth in [t, t + 8] given that 

the system was in state n at time t, 

H n 8 + 0(8) = Probability of a death in [t, t + 8] given that 

the system was in state n at time t. 

Letting Jf (t) be the random variable representing the number of 
circuits at time t, the total probability of a birth in [t, t + 8] may be 
written as 
Prob{a birthing, t + 8]) 

= X Prob{a birth in [t, t + 8] and jV(t) = n) 

= X Prob{a birth in [t, t + 8]\ JT(t) = n}Prob{^(t) = n) 

n 

= £[A„8 + o(S)]P n (*), (61) 

n 

where the certain event and Bayes' rule were successively used. Simi- 
larly, the total probability of a death in [t, t + 8] is £„ [ii n 8 + 0(8)] 
P n (t). With this information, a differential equation for b(t) may be set 
up as follows: 

b(t + 8) = fe(£)[Probability no event or a death] 

+ [b(t) + l][Probability of a birth] + 0(8), 

where o{8) is the contribution of more than one birth. 

b(t + 8) = b(t)[{l - £ (\ n 8 + M + o(8))P n ) 

n 

+ £ ( M „S + 0(8))P n ] 
n 

+ [ b(t) + 1] £ [\ n 8 + o(8)]P n (t) + 0(8). (62) 



• o(-): *?' -» R' is such that lim^„^— = 0. 
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In the limit S — > 0, eq. (62) becomes 



-rb(t)= I \nPn(t). 
at „_n 



(63) 



In a similar manner, the following differential equations for d(t) and 
for m(t) may be obtained: 



Since 



-rd(t)= I tinPM), 
at „=o 

4lMt)= l (Xn-flJPnit). 

at „_o 



m(t) = £ nPn(t), 



(64) 
(65) 



one can use relations (3) to write eq. (65) as a differential equation for 
m(t). Its solutions have already been given in eqs. (4) and (6). With 
knowledge of the mean, eqs. (63) and (64) may be solved for b(t) and 
d(t), respectively. The results are easily found to be 



b(t) = 



A A 



\0t* 



+ (A/io + P)t 



d(t) = 



7 n e(t {e - 1) — 

A A 



\pt* 



where 



and 



+ Xnot 



n e {{ = no + — 
A 



(Af*/i), 

(A = /i), 



(66) 



(67) 



A = X - fi. 

APPENDIX B 

Queue Discipline and Held-Order Probability 

In the model described in this paper, the queue discipline followed 
is Blocked Customers Held (bch). Let the random variable T denote 
the sojourn time of the customer, i.e., the total time he spends in the 
system, either waiting for service or being served. The assumption 
inherent in the bch queue discipline is that the customer will spend 
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time T in the system, after which he will depart, regardless of whether 
he has been served or is still waiting for service. The choice of ju n = 
n (i implies that the sojourn times have a negative exponential distri- 
bution. 

In special services, if the sojourn time distribution is in fact negative 
exponential, then the queue discipline used here should be correct. If, 
on the other hand, it is the service-time distribution that is negative 
exponential, then the bch queue discipline assumed here may still be 
approximately correct if the average waiting time of a customer is 
much smaller than the average service time. 

To compute the held-order probability, care must be given as to 
whether the held order is seen by an outside observer or by an arriving 
customer. For processes with Poisson input, it is well known that the 
distribution P n (t) seen by an outside observer is identical to the 
distribution ir n (t) seen by an arriving customer (see Section 3.2 of Ref. 
3), and hence the distinction is unimportant. For the process described 
in this paper, however, the distributions are different. Define 

P n (t) = Vrob{jV(t) - 7i}, (68) 

TT n (t) — Pmh{Jf(t) = n | a customer arrives at t + ). (69) 

Expression (69) is the probability that a customer who arrives at t 
finds n other customers being served or waiting to be served. Letting 
the event A refer to the arrival of a customer in the interval (t, t + 
8], and using conditional probabilities, one may write ir n (t) as the 
following limit, if it exists: 



1Tn(t) = 



limProb{^U) = n,A) 

6—0 

lim Prob{A} 



= lim 

6—0 



6—0 

Prob{A | Jf(t) = n}¥roh{jV(t) = n) 



£ ¥roh{A\JT{t) = j}Vroh{JT{t) =j) 
/-o 

Since Frob{A \ Jf{t) = n) = \ n 8 + o(5),* one obtains, with the help of 
(68), 

*,(*) - . X " Pn< ' ) • (70) 

S XjPAt) 

y'-o 

The held-order probability may thus be defined as 

h(t) = ^ P n (t) as seen by an outside observer (71) 

n-d+l 



o(-): R 1 -* R 1 »s such that lim5-o^ = 0. 

o 
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or 



Y m n=d XnPn(t) 

h'{t) = -=^ — as seen by an arriving customer, (72) 

Lj= Q x J p At) 

where d is the number of servers. Expression (72) is the conditional 
probability that if a customer were to arrive at t, he would find all the 
servers engaged. This quantity is known in congestion theory as the 
transient call-congestion function. 19 Expression (71) is the probability 
that at least one customer is waiting to be served. This quantity is 
similar but not identical to the transient time-congestion function, 
S(t), which is the probability that all servers are busy at time t, and 
which may be written as 19 

S(t) = i Pn(t). - (73) 

n=d 

For Poisson input (A* = A), it may be shown that 

h'(t) = S(t) > h(t). 
For the Kendall process, the relationship is 

h'(t) > S(t) > h(t). 

For the Poisson input case, since S(t) is equal to h'(t), the time- 
congestion function may be used as a meaningful measure of the held 
orders. For the Kendall process, on the other hand, S(t) does not 
describe the held orders as seen by either an arriving customer or an 
outside observer. Consequently, the time-congestion function is not 
believed to be a meaningful measure of the held orders. Throughout 
this paper, expression (71) was used for the held-order probability, 
although eq. (72) could have been used instead. 

Both h(t) and h'(t) are instantaneous quantities. Since the Kendall 
process with A > fx is not ergodic (i.e., space averaging is different than 
time averaging), a space average must be made when measuring either 
of these quantities. Thus, one cannot measure h(t) or h'{t) by exam- 
ining one sample for a long enough time; rather, one needs an ensemble 
of samples. Because of these measurement difficulties, an open ques- 
tion remains as to whether this type of held-order probability is the 
best measure of the provided service, or whether other quantities such 
as the time average of h(t) or h'(t), or the duration of the held order 
might be more meaningful. Nevertheless, it is clear that eqs. (71) and 
(72) are some measure of the provided service and, as such, are useful 
when comparing different special services provisioning methods meant 
to provide the same level of service. 
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